Modeling the mechano-chemistry of the 029 DNA translocation motor 
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We present a study of the DNA translocation of the bacteriophage </f>29 packaging molecular motor. 
From the experimental available information we present a model system based in an stochastic flash- 
ing potential, which reproduces the experimental observations such as: detailed trajectories, steps 
and substeps, spatial correlation, and velocity. Moreover the model allows the evaluation of power 
and efficiency of this motor. We have found that the maximum power regime does not correspond 
with that of the maximum efficiency. These informations can stimulate further experiments. 



PACS numbers: 87.16.Nn, 87.16. A-, 05.40.-a 

I. INTRODUCTION 

Molecular motors, discovered in living cells almost 30 
years ago, constitute essential ingredients for life, and 
have been extensively studied since then with two ap- 
proaches: biochemical and/or physical Many mo- 
tor are ATP fueled like kincsins and dyneins, which 
can "walk" along a cellular microtubule and perform 
transport features inside cells. Other examples include 
myosins, able to contract the muscles tissues with a 
strongly synchronized movement; or translocator motors, 
which move biopolymers across membranes [5] [5] . Other 
motors are driven by an ion flux such the bacterial nag- 
ella motor (BFM), which propels the bacteria in a fluid 
media. 

Many important cellular processes use the physical 
mechanism of the translocation of biomolecules through 
cell membranes [4j [5] such as proteins, RNA and DNA 
strands. The interest in this subject is twofold: under- 
standing the relevance of this mechanism of living sys- 
tems, but also to investigate the enormous technological 
possibilities of controlling it by using synthetic materials 
able to imitate this important biological function [5] . In 
this sense, the passage of biomolecules through nanopores 
involves many disciplines, from biology to nanothechnol- 
ogy, such as the passage of mRNA through nuclear pores 
[B] or the translocation of DNA in graphene pores [7] . 

In many cases, translocation is driven by constant 
and/or noisy forces inside the pores or by the chemical 
potential difference between both sides of the membrane. 
In other cases, the process is assisted by an intermem- 
brane ATP-based molecular motor [8]. 

To this last category belongs the motor of the 029 bac- 
teriophage which performs the translocation of its DNA 
inside the virus capsid (packaging). This motor translo- 
cates DNA molecules against the huge pressure inside 
the virus capsid, which can achieve 40 atm at the end 
of the process, making it one of the strongest molecular 
machines known. The experimental study of this motor, 
conducted mainly by the Bustamante group [9TU1|. al- 



lowed the knowledge of the DNA dynamics during the 
translocation at a very detailed either time or spatial 
scales. 

Recently a model has been presented in order to de- 
scribe the main features of the 029 bacteriophage translo- 
cation dynamics [12] , as a biological extension of a time 
varying driven translocation of long molecules [131 114) . 
The idea of the model is to drive with a constant force 
a polymer chain in one direction, while in its activated 
state, but it leaves the polymer diffuse freely when the 
motor is inactive [TSJ, [T51 |TB] . The mechanism is the origin 
of the Michaelis-Menten (MM) polymer velocity, related 
to a microscopic re-interpretation of the MM enzymatic 
reaction [12] . 

Although the detailed mechanism of the inner struc- 
ture of the motor remains to be better clarified, a 
sub-stepping feature has been also revealed experimen- 
tally [TT] . These experiments show that each cycle of the 
motor consists of two well differentiated processes. First 
a purely catalytic process in which around four ATPs are 
bound to the pentameric ring. During this phase, the 
DNA does not advance. Hence, this process is referred 
as the dwell phase. The second process of the cycle is the 
burst phase in which the four ATPs are hydrolyzed with 
a total DNA advance of ~10 bp (3.4 nm). Slowing down 
the hydrolysis be means of an external hindering force, 
a finer discretization of the main step is found obtaining 
four substeps of ^2.5 bp (0.85 nm) each one, coinciding 
with the number of ATP hydrolysis processes. 

The present study points its attention to the specific 
sub-step structure of the trajectories, here modeled as 
sub-activation states of the simpler ATP machine de- 
picted in [12l|T6]. As it has been already justified in there, 
for the experimental constraints and the small dimension 
of the motor, the translocation can be approximately de- 
scribed by a rigid chain pushed in a 1-dimensional do- 
main. 

In the next section all the details of the model are de- 
scribed and correlated with the biochemical available in- 
formation. In Section III we present the numerical results 
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with the corresponding analysis and comparison with the 
experiments. Finally we end with some comments and 
conclusions. 



II. THE MODEL 

The full experimental set up consists of the packaging 
motor, the DNA strand and two polystyrene beads of 860 
nm of diameter. However, since the strand is always kept 
stretched, its elastic motion can be neglected. Thus, the 
problem can be considered as a solid motor-bead system 
advancing along a DNA track. The dynamics of such a 
nanoscopic system can be described by an overdamped 
Langevin equation, 



1 x = -V'(x,t)+F E + Z(t), 



(1) 



being 7 the friction coefficient carried in the motion of 
the motor, x the relative position of the motor along the 
track, V(x, t) the motor potential describing the working 
process of the motor, Fe the external force applied on the 
motor and £ (t) the thermal force which can be described 
as white Gaussian noise of zero mean and correlation, 



(2) 



The value of 7 used in the current analysis (Table [T]) 
is an heuristic extrapolation from the observed experi- 
mental data. All the biochemical relevant information of 
the working cycle of the motor, as it is described in Refs. 
[XT] is allocated in the flashing potential V(x,t). Since 
the working of the motor presents two different regimes it 
is expected that two different potentials apply Vd (x) and 
Vb(x) for the dwells and the burst phases respectively. 
Therefore, the temporal and the spatial dependence of 
the potential can be separated as, 

V(x,t) = V D (x) + {V B (x)-V D (x))ri(t), (3) 

where = 0, 1 sets the state of the motor potential. 
r)D = for the dwell phase and 775 = 1 for the burst 
phase. This description introduces a flashing mechanism 
for the state transition and implies that the change in 
conformation occurs in a faster timescale than the actual 
motion of the strand. This difference of timescales has 
been reported for several motors [El [T7] . 

The structure of each potential landscape can be recon- 
structed from experimental results [11 . For the dwell 
state, there must be a minimum for each equilibrium po- 
sition that keeps the DNA fixed when no ATP is present 
in the system and an external force pulls the strand. A 
possible scenario is that of a periodic potential with the 
periodicity of a substep Iq —2.5 bp. The simplest poten- 
tial with such a periodicity is a ratchet potential (Fig. 
[l]). Its height, V^ AX , sets the magnitude of the fluctu- 
ations around the dwell state. However, there are two 



limitations to its value. On one hand the height of the 
potential must introduce a force that is larger or compa- 
rable to the stall force of the motor ~ 80 pN. Conversely, 
a very energetic potential may artificially introduce extra 
energy in the system through the potential flashing. The 
values used in the current work are gathered in Table |l| . 

On the other hand, the excited or burst potential must 
introduce a net force which drives the motor to the next 
minimum of the relaxed potential. This force results from 
the new form of the potential and it can be related with 
a change in the structure of the protein that modifies the 
equilibrium positions of the strand. Since the potential 
is determined by the motor structure, this modification 
is confined to a region of space, i.e. only a region of the 
Vd(x) is modified when switching to Vb{x) (see Figure 
[T]). Note that this scenario is different to the introduction 
of a constant force with no spatial limitations where the 
motion and the energy supplied are uncontrolled. 

The height of the burst potential Vq is not a free pa- 
rameter but is determined by the free energy used from 
the hydrolysis of a molecule of ATP (Fig [lj . Obtaining, 
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Finally, it is necessary to set the temporal dynamics 
of the motor. The easiest description of the dwell time 
to is the necessary time to take irreversibly four ATPs 
independently with a constant reaction rate /catp «-e. 
the stochastic time resulting of adding four exponential 
events. Hence, the resulting probability density is the 
Erlang distribution, 



P(t D ) = 



ATP b D 



3! 



Which has an average time = -j^^ 
randomness 
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and returns a 
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Experimentally, an hyperbolic dependence of the reac- 
tion rate &:atp has been observed an characterized with 
rates &atp and ^atp- This description is in agreement 
with the classic Michaelis-Menten theory and it allows to 
write A: atp as, 



fc A TP 



( ^ATP , .0 

\|atpJ ATp 



(7) 



On the other hand the burst phase can be split in 4 
substeps. Each substep consists of a mechanical process 
(assisting the advance of the strand) plus a catalytic pro- 
cess (that generate the substep dwells). The stroke dura- 
tion of each substep can be computed directly from the 
dissipative dynamics leading the motion, 
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Additionally, another factor remains unknown from 
the trajectory, which is how the substep dwell is extended 
along the two conformational states of the protein i.e. if 
the motor stays in the excited conformation more time 
than t mec h prior to the flash to the relaxed state. This 
time is related to the existence of a process triggering 
the flashing of the potential that has a scale around the 
ms. Such a factor is very important since it will allow 
the motor an additional time to reach the minimum in- 
creasing this way the coupling ratio of the motor. This 
reaction time t r can be introduced as an extension of 
the mechanical time t mec h- For the sake of simplicity, t r 
will be introduced deterministically and the effect of its 
length in the motor performance will be studied in the 
following section. 

On the contrary, the catalytic event after each substep 
has been experimentally characterized as a leading reac- 
tion of a rate k s (FE) that depends on the external pulling 
force applied, and presents a probability density 
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The rate k s at large external pulling forces (-40 pN) is 
measured to be k s = 22 s -1 , while it is very quick at low 
forces and cannot be measured (k s > 1ms -1 ) [TU]. A 
possible description for this rate is an exponential rate 
dependence, 
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with F E ~ -52pN. 

The exact dependence of the rate with the external 
force is not important for the current description as long 
as k s keeps the values mentioned above. With the de- 
pendence in ( 10 ) for low forces, the substep dwell time 
is k~ 1 {F E = 0f = Therefore, i° measures the substep 
dwell time which must be a fast time bellow the millisec- 
ond. Once tg is fixed, the rate that accomplishes that 
k s (— 40pNnm) = 22s -1 can be calculated. 

This closes the drawing of the dynamics of a cycle (Fig. 
[T]), which it is worth to summarize here. It starts in the 
dwell state rj(t) = and switches to the burst state after 
a stochastic time to- After that, the potential starts the 
first stroke of the cycle flashing to the burst state rj(t) = 1 
that lasts a deterministic time t mec h + t r flashing back to 
the dwell potential for a stochastic time t s . This process 
is repeated three more times finishing the burst phase 
and therefore the full cycle of the motor. 

The value of the parameters used has been summarized 
in Table I. 



Figure 1: (Color Online) Flashing potential describing the 
motor cycle. The cycle starts with the dwell phase (orange). 
Once the four ATP are bound, the burst phase (blue) begins. 
The cycle is not fully depicted, in order to complete the cycles 
two additional substeps are necessary. 



Parameter 


Value 


Jo 


0.85 nm 


7 


200 pN ms/nm 


kT 


4.1 pN nm 


AGatp 

i/MAX 
V D 


90 pN nm 
150 pN nm 


K ATP 


650 ms (M 


I AXP 


28.8 ms 


F s 


6pN 




0.05 ms 




1 ms 



Table I: Parameters used in the simulation. 



III. RESULTS 

Introducing the dynamics and the potential landscape 
in eq. ([!]) and using the parameters of (Table [TJ) the 
description is complete and the stochastic equation of 
motion (JlJ can be simulated. 

The resulting trajectories are steplike with a sub- 
stepping dynamics dependent of the external force. In 
Fig. [2] it is possible to notice that when the pull force is 
low {Fe = — 8pN), then the substep feature of the tra- 
jectory is hardly visible, while conversely it is well visible 
for pulling force of Fe — — 40pN. The statistical proper- 
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Figure 3: Spatial correlation of the two trajectories in Fig. 
[5] for a run of 10 seconds. F E = — 8pN (top) and — 40pN 
(bottom). 

ties of the substepping can be better studied through a 
correlation analysis of the trajectories (Fig. [3]), obtained 
as the autocorrelation of the spatial density profile of the 
motor P(x), 

C(Ax) = (P(x + Ax)P(x)) x , (11) 

where P{x) is computed as the position histogram for the 
trajectories. The correlation reveals the spatial stepping 
behavior described. Again, the substepping features are 
very clear for the strongest pull force, and quite absent 
for the lowest value. The evidence of the stepping mecha- 
nism for higher pull forces coincides with the experimen- 
tal observations. The two trajectories in Fig. [2] can be 
compared with the experimental results in Ref.[20 , Figs, 
lb and 3a for Fe = — 8pN and Fe — — 40pN respectively. 
Analogously, the two lines of Fig. [3] represent quite finely 
the plots in Fig. lc and 3b of the same article. 

The trajectories also show that the greater the external 




Ax(bp) 

Figure 4: Spatial correlations for three values of t r , namely 
t r — 0, 0.1, and 1ms. The pull force applied is Fe = — 8pN. 



opposing force is, the slower is the motor. There are two 
reasons for this effect. On one hand a greater hindering 
force reduces the net stroke force of the motor. On the 
other hand it tilts the potential increasing the number 
of failed steps. These are the steps where the full ATP 
hydrolysis energy is consumed but the motor does not 
arrive at the following minimum of Vd(x). This effect 
can also be studied through the deterministic time t r 
observing that an increase in the reaction time prior to 
the flashing to the relaxed state brings an increasing of 
the coupling ratio (Fig. [4]). It is interesting to note that 
the time t r in all the cases studied here is very small and 
practically does not increase in any case the time of the 
cycle. 

From the trajectory, the dependence of the average 
velocity of the motor v with the external force F E can 
be studied obtaining the classical decaying force-velocity 
curve (Fig. [5]). The same occurs when the dependence of 
the velocity with the ATP concentration is studied (Fig. 
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Figure 5: (Color Online) Average velocity of the strand ver- 
sus the external hindering force for different concentrations of 
ATP. 
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the motor and ^ AVa as h are all the potential flashing in- 
crements corresponding to the same number of hydrolized 
ATP's of this trajectory [IB]. 

In Fig. [8] one can see the characteristic curve of the 
efficiency with its maximum at F E ~ — 52pN, a value 
close to the stall force. The maximum of efficiency is ~ 
38%, less than most of ideal motors where the efficiency 
is 50%. It is worth to remark that the efficiency does 
not depend very much on the ATP concentration and 
that the maximum of efficiency does not coincide with 
location of the maximum power. 
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Figure 7: (Color Online) Power of the motor as a function 
of the external force for different concentrations of [ATP]. 



IV. COMMENTS AND CONCLUSIONS 



Figure 6: (Color Online) Average velocity of the strand ver- 
sus [ATP] for different values of the external force. 



[6]) where Michaelis-like curves are seen. The latter evi- 
dence has been also found in Ref . [I2j , and experimentally 
reported in Ref. [TO] , 

The analysis of the dynamics allows also to study the 
energetics of the model: power and efficiency, where the 
the power extracted form the motor is defined as F E (v) . 
Fig. [7] shows the classical parabolic shape of the power 
that becomes null in two extreme points, at F E = 0, and 
at the stall force of the motor. This behavior points out 
an optimum force for which the power extracted from 
the motor is maximum. It is interesting to note that the 
curve is not symmetric. Furthermore, the external force 
needed for a maximum power changes with [ATP] . 

The efficiency of this motor can also be calculated as 
the ratio between the output work and the input energy 
through ATP consumption, 
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In this work we have introduced a model which can 
explain the trajectory of the 4>29 packaging motor, and 
in particular its substep features, using an on-off sys- 
tem (flashing potential) which mimics the ATP hydroly- 
sis processes to produce mechanical work. 
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being AL the total distance advanced in a trajectory of 



Figure 8: (Color Online) Efficiency of the motor versus the 
force Fe and two ATP concentrations 
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The model is based in a careful analysis of the chemical 
process that occurs in the inner mechanism of the motor, 
evidenced by the known structure reported in Ref. [TT] . 
The inner mechanism has been modeled here as four ex- 
ponentially distributed ATP waiting times, followed by 
the subsequent ATP hydrolysis events, which give rise to 
the detailed trajectory sunsteps. 

The resulting spatial correlation calculated with dif- 
ferent pull forces applied to the chain, shows the typical 
trend revealed experimentally. Also the velocity of the 
translocation process for different ATP concentrations, 
shows a clear agreement with the experimental outcomes, 
confirming the Michaelis-Menten dependence as a result 
of the ATP machine functioning. 

The authors of [TT] report as surprising the fact that 
the correlation distribution decays with Ax. In our 
model this result appears by itself and its origin is due 
to missing steps which induces a lost of the correlation 
for long distances. 

Another interesting point to remark is the presence 
of high fluctuations observed in the experiments if com- 
pared with the ones of the simulations. A possible expla- 
nation is that we simulate a rigid system formed by the 
motor, the DNA chain and the beads with an effective 
friction constant. Conversely the chain is not rigid, and 
this may soften the transmission of the motor dynamics 
to the bead. 



As in other investigations [T2TU4| . we concentrate our 
analysis on a one-dimensional model. This approach is 
justified because the polymer translocation in the cited 
experiments is usually performed with optical traps, and 
so the polymer is maintained stretched during the pro- 
cess. The polymer is considered rigid, in order to check 
the features of the machine more than the dissipation 
given by a smooth chain. 

That way we have phenomenologically described the 
details of the bacteriophage trajectories, using assump- 
tion on the motor functioning only. Actually, the model 
can be extended to other akin motors, such as the re- 
cently the study of the packaging motor of bacteriophage 
T4 |21j . For this motor, pause-unpacking steps in low 
ATP concentrations and large forces were experimentally 
observed . Our model can address this issue with the ap- 
propriate changes in the biochemical parameters. 
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